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A SIMPLE MODEL FOR THE INTERACTION OF PLANETARY 
AND SYNOPTIC-SCALE WAVES IN THE TROPOSPHERE 

by 

Lee E . Branscome 

Associate Professor of Meteorology 
and Physical Oceanography 
University of Miami 
Miami, Florida 

ABSTRACT 


The time-mean circulation of the troposphere during the 
Northern Hemisphere winter is dominated by planetary- 
scale stationary waves superposed on longitudinal -mean 
westerly flow. The movement and generation of weather- 
making, synoptic-scale transient waves are controlled by 
this large-scale circulation. The interannual variability 
of planetary-scale waves determines the severity of the 
winter season over broad geographic regions through such 
factors as frequency of storms, precipitation and tempera- 
ture. 

The transient waves or eddies play an important role in 
maintaining the heat and angular momentum balance of the 
atmosphere. By maintaining the westerly flow in mid- 
latitudes, these waves are responsible for the mean zonal 
flow over large-scale topography that forces, in part, the 
stationary planetary waves. In addition, the precipitation 
and subsequent latent heat release in the organized "storm 
tracks", i. e. , preferred regions of transient wave activity, 
form a longitudinally-varying pattern of diabatic heating 
that also forces stationary waves. On the other hand, the 
transient, baroclinic waves damp the planetary-scale waves 
by the down-gradient flux of sensible heat. 

A conceptual model that incorporates these feedbacks 
between the two scales is developed and guides the formu- 
lation of a quantitative low-order model which retains a 
few basic physical processes and wave components. In the 
low-order model the fluxes of sensible heat and momentum 
and latent heat release by transient eddies are parameter- 
ized in terms of the evolving large-scale circulation. 
Radiative forcing and damping, planetary-scale variations 
in moisture and topography and Ekman dissipation are also 
included. Testing of this model is in progress. The results 
will help interpret experiments from sophisticated gene- 
ral-circulation models and ascertain the need for satellite 
observations of precipitation. 
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INTRODUCTION 


The time-mean circulation of the mid-latitude tropos- 
phere during the Northern Hemisphere winter is dominated by 
westerly (eastward) flow and stationary waves of planetary 
scale (see Fig. la) . The movement and generation of synop- 
tic-scale (1000-4000 km) , high-frequency (2-6 day period) 
transient waves are controlled by the large-scale circula- 
tion (see Fig. lb and Background section below).' These 
synoptic-scale waves, which have significant circulation 
throughout the troposphere, are manifested at the surface 
by the mid-latitude storms or cyclones that determine the 
daily weather. 

The planetary-scale waves have substantial interannual 
variability and their positions and amplitudes averaged 
over a given winter can deviate significantly from the 
climatological norm (see Fig. 2) . Since these waves direct 
weather-making cyclone activity, the wintertime precipita- 
tion and temperature distributions can also deviate from 
normal patterns, thus determining the severity of the win- 
ter season over broad geographic regions. 

Therefore, it is of considerable importance to under- 
stand the mechanisms responsible for forcing the planetary 
waves to improve weather prediction, seasonal forecasts and 
climate models. Because the transient waves have transport 
properties that modify the global-scale circulation, we 
suspect that they play a crucial role in the dynamics of 
the planetary waves. We anticipate that feedbacks or 
interactions between the two scales of motion continuously 
operate . 

One method for evaluating the controlling feedbacks 
relies on simple "low-order" models which retain only a few 
physical processes and wave components. These models allow 
us to isolate the dominant mechanisms that control plane- 
tary wave behavior. This approach to the problem will be 
developed in this study. 
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OBJECTIVES 

The objectives of this study were to: 

1) develop a conceptual model for the generation of plane- 
tary-scale waves and their interaction with synoptic-scale 
waves or eddies 

2) develop an idealized, quantitative low-order model to 
simulate the basic feedbacks of this system 

3) develop a computer code to solve the equations of the 
quantitative model 

4) compare the modeling with observed properties of the 
atmosphere and ascertain the need for space— based observa- 
tions of certain atmospheric parameters related to these 
feedbacks or interactions. 











ORIGINAL PAGE IS 
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Fig. 2. a) Mean geopotential height for January 1981 at 500 
mb. Note amplified ridge-trough-ridge pattern over N. 
America-Atlantic . b) As in a. except for 1000 mb (similar 
to surface pressure). Arrows show direction of near-surface 
flow over N. America. This pattern led to unusually cold 
January in the Eastern U.S. with warm, dry weather along 
the Pacific coast of N. America. (from Wallace and Blackmon 
1983). 
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BACKGROUND 


1. The role of transient eddies in the global circulation 

The latitudinal variation of incoming solar radiation, 
primarily in the visible wavelengths, produces a latitudi- 
nal temperature gradient with the tropics being warmer than 
high latitudes. This contrast in heating introduces 
"available" potential energy into the atmospheric system 
and provides the energy for atmospheric and, consequently 
also, oceanic motions. As with many fluid systems with 
spatial contrasts of some parameter related to energy 
(here, the relationship is between temperature contrast and 
the potential energy in a hydrostatic fluid) , eddy motions 
appear, esp. through an instability process, and attempt to 
eliminate the gradient of the energetic parameter. 


In the case of the atmosphere, this instability process 
is called baroclinic instability since it relies on the 
baroclinic nature of the fluid. The differential solar 
heating enhances the baroclinicity (essentially, the tem- 
perature contrast along a constant pressure surface) , while 
the baroclinic eddies or waves attempt to eliminate the 
thermal contrast by transporting sensible heat down the 
temperature gradient. The waves derive thei f irciu 

the energy which is introduced into the system by the 
insolation. 


The overall heat balance of the global atmosphere is 
depicted in Fig. 3. Relatively large amounts of solar 
radiation are absorbed by the earth in low latitudes. The 
baroclinic waves, which are the same transient waves dis- 
cussed in the Introduction, transport heat poleward. At 
high latitudes the energy is emitted to space in the form 
of infrared radiation. These high latitudes are actually 
warmer than they would be if no heat was transported pole- 
ward and the temperature was determined only as a result of 
a purely radiative balance at those latitudes. In fact, 
these latitudes emit more radiation (in the form of 
infrared radiation) than they absorb (in the form of solar 
radiation) . 

An additional characteristic of the baroclinic insta- 
bility process is that these waves also transport angular 
momentum poleward. This transport plays an essential role 
in the maintenance of the observed latitudinal distribu- 
tion of surface winds. Fig. 4 shows the basic angular 
momentum balance of the earth-atmosphere system. Middle 
(low) latitudes have westerly (easterly) surface winds. 
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Incoming Solar Radiation 


Fig. 3. Schematic diagram of global heat balance 


Pole 



Angular Momentum 
Transport 


Fig. 4 Schematic diagram of atmosphere-Earth angular 
momentum balance 


Because the atmosphere is rotating faster (slower) than the 
earth in middle (low) latitudes, it loses (gains) angular 
momentum to (from) the earth through frictional drag at the 
surface. The bar aclinic waves are partly responsible for 
this surface wind distribution by providing the primary 
mechanism for the transport of angular momentum into mid- 
latitudes. 

In addition, baroclinic eddies also transport water 
vapor poleward. From experience we know that southerly 
winds in advance of a cold front and associated cyclone are 
warm and humid, whereas the northerly winds behind the 
front are cold and dry. Averaged over the entire weather 
system, this correlation implies the northward transport of 
water vapor. The source of the water vapor is primarily 
evaporation from the subtropical ocean surface. By the 
Clausius-Clapeyron relation, the saturation vapor pressure 
increases rapidly with air temperature, allowing the warm 
southerly flow which obtains moisture from the subtropics 
to carry large amounts of water vapor poleward. 

Part of the baroclinic instability process includes a 
correlation between upward motion and relatively warm air. 
This correlation allows potential energy existing in the 
thermal contrasts of the wave field to be converted into 
eddy or wave kinetic energy. A consequence of this uplift- 
ing of warm air in a moist atmosphere is adiabatic expan- 
sion and cooling of air parcels and, inevitably, condensa- 
tion and precipitation. This precipitation is most evident 
in the cyclonic (counterclockwise rotation in the Northern 
Hemisphere) portion of the synoptic-scale wave. This con- 
densation involves the release of latent heat and, there- 
fore, introduces an additional heating mechanism, one that 
is dependent on the activity of the transient baroclinic 
waves . 

We saw in Fig. lb that the synoptic-scale, transient 
wave activity is largely confined to "storm tracks" over 
the northern oceans. Very reliable estimates of precipita- 
tion over the oceans are not available. We expect, how- 
ever, that precipitation and latent heat release are maxim- 
ized in these oceanic storm tracks, esp. since a ready 
source of moisture from the subtropical ocean is present. 
Thus, the maxima in latent heat release in these storm 
tracks introduces longitudinal asymmetries in diabatic 
heating. 


In summary, baroclinic wave activity controls 
the strength and sign of the surface wind in a long-term 
mean. In addition, baroclinic waves extract potential 
energy from the larger-scale circulation and, in doing so, 
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act to eliminate horizontal temperature contrasts of the 
larger scales. The storm systems of these waves produce 
precipitation and planetary-scale variations in diabatic 
heating averaged over several storm occurrences. 

2. "Storm tracks" 

As seen briefly in the Introduction and Fig. lb, the 
planetary-scale circulation controls the movement and gen- 
eration of the synoptic-scale waves. The nature of this 
control is not well -understood theoretically, but some gen- 
eral comments from observations and basic theory are pos- 
sible. 

First, the behavior of mid-latitude atmospheric waves 
with horizontal wavelengths of about 1000 km or longer is 
governed by "quasi-qeostrophic" dynamics (see Appendix for 
brief description of the quasi-geostrophic system) . The 
phase speed of a linear wave in the presence of a back- 
ground mean zonal (east-west) flow can be derived from the 
conservative forms of (A. 1-3): 


c, phase speed = u - ( 1 ) 

k 2 + l 2 + m 2 f 2 /N 2 


where u is the mean zonal velocity, k,l,m the zonal, lati- 
tudinal and vertical wavenumbers, respectively, and /3 
the latitudinal derivative of the Coriolis parameter which, 
along with the remaining variable, N, is defined in the 
Appendix. 

From (1) we see that as the wavenumbers increase, the 
wave phase speed is closer to the mean zonal velocity since 
the second term in (1) , which is always positive, becomes 
smaller. The shorter waves are more likely to be "steered" 
by the large-scale velocity field. Therefore, as evidenced 
by Fig. 1, the synoptic-scale waves are guided in their 
movement by the planetary-scale circulation. The planetar- 
y-scale waves are more likely to be slowly propagating or 
stationary. The second term which involves the wave's 
advection of planetary vorticity is larger for these waves 
and tends to cancel the mean flow advection of the wave 
vorticity as represented by the first term in (1) . 

An additional consideration is the generation of synop- 
tic-scale disturbances. As mentioned above, they obtain 
their energy from the large-scale temperature contrasts. 
These contrasts are maximized in long-term time means at 
the base of the planetary-scale troughs. Thus, transient 
synoptic-scale disturbances appear to be preferentially 
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generated in these regions and are then carried downstream 
by the large-scale velocity field. 

3 , Mid-latitude variability 

Some observational and modelling analyses (e.g., Horel 
and Wallace, 1981; Pitcher et al . . 1983) suggest that lar- 
ge-scale, mid-latitude circulation anomalies are strongly 
influenced by changes in atmospheric forcing from the 
earth’s surface, esp. through fluctuations in tropical sea 
surface temperatures (SST's). SST anomalies from the cli- 
matological norm, esp., the elevated temperatures during 
the El Nino event, alter the pattern of tropical atmo- 
spheric heating and may induce considerable changes in the 
mid-latitude large-scale circulation. 

This hypothesis omits the significant interaction of 
synoptic and planetary-scale waves in mid-latitudes. 

Recent modelling and observational studies indicate that 
persistent circulation anomalies, which last for two weeks 
or longer and are fixed in longitude (so-called "blocks"), 
are initiated, maintained and, perhaps, destroyed by tran- 
sient, synoptic-scale disturbances (Dole, 1983; Mullen, 
1985; Colucci, 1985). Indeed, the large-scale circulation 
often shifts dramatically from one persistent pattern to 
another within a few days. This rapid transition cannot be 
easily explained by the slow evolution of tropical SST's 
which transpires over a few months or more. 

4. Directions for modelling interaction of synoptic and 
planetary waves 

It will be our premise that synoptic-scale waves conti- 
nuously interact with the planetary-scale circulation. 

To successfully model the planetary-scale waves over a few 
days or on a climatic basis, the effects of the synoptic- 
scale eddies must be included. In the following section we 
will consider mechanisms which are responsible for the gen- 
eration of planetary-scale waves. Accounting for the 
nature of synoptic eddy transports as discussed above, we 
will develop a conceptual model for the interaction of 
these two scales of motion. 

Our direction in this project will be to understand the 
interactions of the two scales in a climatological sense. 
That is, we are interested not in a short-range prediction 
of the circulation, but, instead, in what mechanisms oper- 
ate in the time mean to maintain the observed pattern of 
stationary planetary waves seen in Fig. la and the time- 
mean transports by the synoptic-scale eddies. 
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CONCEPTUAL MODEL 


1. Mechanisms for generating planetary waves 
a . Topography 

Two mechanisms for generating planetary-scale waves are 
considered here in the context of the preceding discussion. 
First, steady surface flow over large-scale topography 
(e.g., Rockies, Himalayas) will generate stationary plane- 
tary-scale waves. This concept has been treated exten- 
sively (one of the earliest papers was by Charney and 
Eliassen, 1949; a recent excellent review is given by Held, 
1983) though not in the context of interaction with synop- 
tic-scale eddies. 

Fig. 5 shows the basic idea in terms of the quasi- 
geostrophic vorticity equation (A.l) which accounts for 
changes in the absolute vorticity of a fluid parcel in the 
presence of vertical motion. Westerly flow over a mountain 
barrier decreases (increases) the vorticity on the upwind 
(lee) side by the vertical contraction (stretching) of ver- 
tically-oriented vortex tubes. The reduced (increased) 
vorticity on the upwind (lee) side creates an anticyclonic 
ridge (cyclonic trough) , forming a standing wave pattern. 

The forcing for this wave depends on the strength of 
the vertical velocity created by the surface flow over the 
topography. Specifically, the vertical velocity (i.e., the 
velocity component parallel to the gravitational force) at 
the earth's surface is given by the kinematic boundary con- 
dition: 


W(z=h) = V h * 7 h (2) 

where h=h(x,y) is the topographic height and Vjj is the 
surface horizontal velocity. - 

As noted in the Background section, the time-mean sur- 
face wind distribution is largely controlled by the angular 
momentum transport by baroclinic transient (i.e., synoptic- 
scale) waves. Therefore, these waves have an indirect 
influence on the strength of the topographically-forced 
stationary waves. Thus, a key element of any interactive 
model will be this angular momentum transport. 

b. Diabatic heating 

Longitudinal asymmetries in diabatic heating, as repre- 
sented by the non-conservative term, Q, in (A. 2), can also 
induce stationary planetary waves. In the presence of a 


IX-10 



MECHANISMS FOR GENERATING PLANETARY WAVES 


Flow over Topography 



Fig. 5. Forcing of standing planetary wave by steady 
westerly flow over large-scale topography 


wave-like perturbation of Q, the dominant balance in 
the thermodynamic equation (A. 2) for time-mean planetary 
waves is between mean zonal advection of the wave's temper- 
ature pertubation and the diabatic heating (Held, 1983) . 
This balance is depicted in Fig. 6 which shows that the 
temperature field leads the diabatic heating field by 
roughly 90°. 

The thermal-wind equation (A. 3) for this hydrostatic, 
geostrophic motion shows that a relatively cold region is 
associated with a decrease of the streamfunction with 
height. For a wave-like perturbation this implies an 
increase of the Laplacian of the streamfunction (i.e., the 
relative vorticity) with height. Thus, a trough or cyclonic 
circulation intensifies with altitude in a cold region. 

The appearance of longitudinal asymmetries in precipi- 
tation in mid-latitudes encourages the formation of sta- 
tionary planetary waves through this mechanism. The 
strength of this forcing is dependent on the longitudinal 
distribution of moisture and cyclone (synoptic-scale wave) 
activity. Again, interaction between the two types of 
waves is possible through a indirect means. 

c. Damping by synoptic-scale eddies 

As suggested above, synoptic-scale eddies play a cru- 
cial role in the generation of planetary waves. On the 
other hand, these waves act to destroy temperature con- 
trasts and will damp the planetary waves through the heat 
flux occurring in the baroclinic instability process. 

This effect must also be incorporated in an interactive 
model. The observational analysis of Holopainen (1983) 
indicates that this damping mechanism is an important means 
of limiting the planetary-wave amplitude. 

2 . Feedback diagram 

Fig. 7 shows a schematic diagram which describes the 
means by which planetary waves are forced and interact with 
the transient baroclinic waves. It is not intended to be a 
complete diagram of the general circulation processes of 
the atmosphere, but many important elements are included. 

To a limited extent, one can picture the feedbacks as 
exchanges of energy from one form or one horizontal scale 
to another. The arrow directly connecting the stationary 
planetary waves with synoptic transient waves represents 
the loss of energy to these smaller scales through the 
thermal damping of the baroclinic instability process. 
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MECHANISMS FOR GENERATING PLANETARY WAVES 
Diabatic Heating 


Dominant balance in thermodynamic equation for these waves: 
[u] /cr9' \ = Q 

l e 0 / 

advection of temp, wave heating/cooling 

by westerly flow 



Cold regions (i . e. , 0 1 < 0) ===> ===> / >0 

0 O 5 z bz 

Relative vorticity increases with altitude (i.e., cyclonic 
circulation intensifies with height) . 


Fig. 6. Forcing of planetary waves by longitudinal asymmet- 
ries in diabatic heating 
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FEEDBACK DIAGRAM 


Differential 
Solar Heating 


Eq.-to- Pole Temp, 
Gradient 


poleward 


Transient Waves 

heat flux 


(baroclinic instability) 


Polevard Angular 
Momentum Transport 


Surface Westerly Plow 
Over Topography 


Planetary Waves 


Storm tracks 


Frictional Dissip 


Frictional dissip,/ 
Radiative damping 


Longitudinally-varying 

y 

Longitudinal var. 

Diabatic Heating (precip.) 

— 

of moisture source 


Fig. 7. Feedback diagram summarizing conceptual model 
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QUANTITATIVE » LOW-ORDER” MODEL 


1. Related studies 

The approach will follow similar low-order models which 
retain only a limited number of horizontal and vertical 
modes in a spectral representation of the evolving fields. 
For example, Charney and Straus (1980) used a low-order 
model to examine the nonlinear equilibrium states of the 
zonally-averaged zonal flow and a single zonal wave that 
exist in presence of radiative forcing, Ekman dissipation 
and planetary-scale topography. An extension of their study 
(Reinhold and Pierrehumbert, 1982) included an additional, 
shorter wave and showed that synoptic-scale waves play an 
essential role in the types of planetary-scale wave config- 
urations that can occur. 

The formulation here is similar to a study by White and 
Green (1982) who also examined the forcing of planetary- 
scale waves in the presence of transient eddy fluxes. Like 
this project they also parameterized the effects of the 
transient eddies rather than explicit ing calculating 
individual synoptic-wave events. However, serious problems 
arise as a result of their physical modelling. 

They exclude diabatic heating caused by latent heat 
release in transient eddies and specify a heating wave 
which is independent of evolving flow parameters (i.e., 
non-interactive forcing) . We anticipate considerable 
interaction between the planetary-scale circulation and the 
transient eddies responsible for the latent heating. 

Their parameterization of the eddy sensible heat flux 
is quasi-linear and diffusive and damps the shortest modes 
most strongly. This mechanism reduces the amplitudes of 
waves which are primarily forced by topography and allows a 
dominance of the longest modes (planetary wavenumber one) 
which are primarily forced by diabatic heating. An improved 
parameterization is used here which depends on the insta- 
bility characteristics of the large-scale circulation and 
does not preferentially damp shorter scales. 

Lastly, radiative heating is specified in their model 
and is independent of the temperature field. Therefore, no 
radiative damping of the planetary-scale thermal waves is 
included. In addition, poleward dynamical heat flux is 
effectively predetermined by the specification of radiative 
heating, prohibiting interaction between transient eddy 
heat flux and the radiative heating. 
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2 . Basic equation and spectral representation 

We will combine many aspects of the conceptual model in 
a quantitative "low-order" model of the general circula- 
tion. First, we use (A. 1-3) to find a single equation for 
r , i.e., the quasi-geostrophic potential vorticity equa- 
tion, 

d { 7*^+ f + f 2 5 ( N ” 2 ^)? = + f Jb (QN“ 2 ) (3) 

dt *5z ' Iz 1 * dz 


The terms on the right-hand side of (3) represents 
sources/sinks. Note that, in the absence of these terms, 
the quantity in the brackets (quasi-geostrophic potential 
vorticity) is conserved with the motion. 

Eq. (3) is a nonlinear PDE and has no general closed- 
form, analytic solutions. First, we model the vertical 
variation of ^ in (3) by assuming 


^(x,y,z,t) = ^ B (x,y,t) F 0 (z) + ^ T (x,y,t) Fi(z) 


where the functions Fq and F]_ are shown in Fig. 

z 



z=H 


IF 0 


z = 0 
F(Z) 


8 . 


(4) 


Fig. 8 Vertical modes, F^(z) 

The barotropic component of the streamfunction, • \ B , has no 
temperature gradients associated with it, whereas the bar- 
oclinic component, rilf h as horizontal temperature gra- 
dients, even at the upper and lower boundaries of the 
model. This latter feature makes this choice of vertical 
modes more realistic for atmospheric modelling than the 
standard two-mode approach (used primarily in ocean models) 
which has no boundary temperature gradients . This modal 
representation produces two equations for and that 
are mathematically analogous to the standard two— level 
model of Phillips (1954) that is widely used in atmospheric 
modelling. The modal formulation has the advantage of 
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easily projecting vertically continuous parameterizations 
of transient eddy heat and momentum fluxes and latent heat- 
ing onto each vertical mode. It also ensures that energy is 
conserved in the presence of topographic forcing and is 
lost through ttkman damping (properties not necessarily 
guaranteed in a two-level model) . 


By substituting (4) into (3) and projecting the equa- 
tion onto the two vertical modes, we find evolution equa- 
tions for y-Q and Upper and lower boundary terms, 

which appear upon taking the projections, can be expressed 
in terms of the vertical velocity at the boundaries. We 
assume that W(z=H)=0 and W(z=0) will be determined by flow 
over topography and Ekman pumping. These equations are 


a 

V X ^B + 

J (^ ^B) + J (^T' 

7 + 

B*lv, 

S5 

at 




“ax 




< & > - fwro) 

H 




•a 

(7 l / T - 

+ J (^r rV ^ 3 ) 

+ J (^B' 



at 




«■ t 

p ax 

= - fW(0)F!(0) + 

< F^> 

- (f/N) 2 <QdJ? 1 /dz> 


(5) 


( 6 ) 


where J is the x-y Jacobian, < > indicates vertical average 
and A* =12.4 (f/NH) 2 . 

The horizontal structures of ^3 and are repre- 
sented by an eigenmode expansion in a periodic channel 
centered in mid-latitudes: 

SP i (t) Gi(x,y) ©i(t)G-[(x,y) 

where the eigenmodes are found from the equation 
V 4 G^ = -a ^ 2 G^ with G^ periodic in x over the length L*, 
dG^/dx=0 and d[Gj[]/dy =0 at y=0,Ly. The brackets indicate 
a zonal (x) average. The latter two conditions ensure no 
meridional velocity or stress at the northern and southern 
sidewalls of the channel. 


We truncate the series in G to resolve only the longi- 
tudinal mean flow and planetary scales. The effects of 
synoptic-scale waves will be parameterized. For initial 
calculations we retain only 6 of the gravest orthonormal 
modes : 

Gi = 2 cos y* 


G 2 = 2 cos 2 y' 



G 5 = 2 sin nx cos 2 y' Gg = 2 cos nx sin 2 y' 

where y' = /Ty/Ly and n-2/V/lj^. In effect, this truncation/ 
leads to a 12 -parameter model for the six components of r -q 
and 

3 . Vertical velocity at the lower boundary 

The vertical velocity, W(z=0) , is given by (2) , 
representing the topographic forcing, plus an additional 
"Ekman pumping" term. The latter term is due to vertical 
motions induced by convergence in the frictional boundary 
layer near the surface. Convergence (divergence) in the 
boundary layer occurs when the interior geostrophic flow 
has cyclonic (anticyclonic) vorticity. The convergence 
(divergence) induces upward (downward) motion at the top of 
the boundary layer which decreases (increases) the interior 
geostrophic vorticity by the vertical contraction (stretch- 
ing) of vortex tubes in the interior (see (A.l)). The 
total vertical velocity is , 

W(z= 0 ) = J( z=0), h) + k 7 r(z**o) 

where k=(A v /2f) 1 / 2 and A v is the eddy viscosity. 

4 . Radiative heating 


We parameterize radiative processes by linear relax- 
ation to the temperature of radiative equilibrium, T e , 
i.e., the temperature field in the absence of dynamical 
heat transports. T e is assumed to be independent of longi- 
tude (x) with its largest value at the southern boundary, 
decreasing montonically to the north. The difference in T e 
across the channel is specified and determines the strength 
of the differential solar heating. Larger values will 
force stronger dynamical heat transports. However, as the 
transports increase, the temperature contrast across the 
channel will decrease and reduce the energy available to 
the baroclinic transient eddies, hindering further poleward 
heat flux. The increased temperature in high latitudes 
leads to stronger radiative cooling and a restoration of 
the meridional temperature gradient which drives more eddy 
heat transport. Thus, the heat flux and radiative pro- 
cesses form an interactive system. 

The zonal-mean temperature field (equivalently, the *P ^ 
field; see (A. 3)) relaxes to T e in the absence of dynamical 
heat fluxes on a time scale of 2z , usu. taken as 20 days. 
The thermal field of the planetary wave is radiatively 
damped on the same time scale. These radiative terms 
appear in the source term, Q. 
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5 . Topography 

As a initial test, we limit the topography to the G 4 
mode. That is, we take h(x,y) =h 0 G 4 (x,y) where 4h 0 is the 
height difference between the lowest and highest points, 
typically, 3-4 km. Later versions of the model will 
examine topographic forcing on various scales. 

6. Parameterization of transient eddy fluxes 

a. Eddy source terms 

The remaining contributions to the source/ sink terms in 
(5) and (6) originate from unresolved transient-eddy acti- 
vity. Those terms are 

= - V* (v'j* ) 


Q = “V* (v'i^'/3z) + Q c 

where $ is relative vorticity and ( j is a time average 
over a period long compared to the synoptic time scale (2-6 
days) and the prime is the deviation from that average. 
Thus, the correlations represent synoptic-eddy fluxes of 
vorticity and heat. The quantity, Q c ,is the latent heating 
due to eddy activity. 

b. Heat flux 

The approach here will follow Branscome's (1983) para- 
meterization of meridional eddy heat flux which is based on 
baroclinic instability theory on a mid-latitude beta-plane. 
This parameterization is generalized to a two-dimensional 
heat-flux vector. Observing the down-gradient nature of 
the heat flux, we write this formula as 

v 1 0 1 = -{AgNH 2 /(0f 2 ) } I V0 1 Ve (1+Y)" 2 exp(-z/d) 

where Y = (?H0N 2 / (fg|V0 1 ) , a dimensionless related to the 
instability of the baroclinic flow; d= H/(l+Y), the depth 
of heat flux; and A=0.6, a correlation constant determined 
by instability theory and/or observations (see Branscome, 
1983) . The parameters inside the brackets are assumed con- 
stant in the model, while the remaining factors vary in 
x, y, and t as the thermal gradient (i.e., p T ) in the , 
planetary-scale circulation evolves. After relating 0 to r T 
this parameterization enters (6) through the Q term. 

Though the heat flux is down-gradient, this formulation is 
highly nonlinear (not linear diffusion) and projections of 
the full flux field onto each G^ mode of ^ T must be made. 

c. Vorticity flux 

The primary influence of the vorticity flux in the 
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dynamics of the planetary-scale circulation comes through 
the meridional transport. This component is responsible 
for the poleward angular momentum transport which drives 
the surface zonal wind. The zonal flux acts to damp the 
planetary-scale stationary waves, but it is less effective 
than the transient eddy heat flux (apparent in analysis of 
Holopainen, 1983) . Therefore, we model only the zonal-mean 
meridional transport, i.e., we take 

^ = - 3[v« $ ' ]/»y. 

It can be shown from the definition of the potential 
vorticity that 

(v*q ] = [v*S*l + (f/N ) 2 a[v* 3 ^*/»z]/*z (7) 

where * indicates deviation from the zonal mean and q is 
the potential vorticity as defined in ( 3 ). Furthermore, we 
find [v 5 ] = - 3[u v*]/3y. Integrating (7) in y and z and 
imposing the condition that [u v ]=0 at y=0,Ly, we find 

I J Ly [v*q*] dy dz = f** (f/N) 2 [v*a^*/3z] | H dy. ( 8 ) 

0 0 * 0 0 

That is, the volume integral of potential vorticity flux is 
proportional to the difference of the upper and lower 
boundary integrals of the heat flux. 

Next, we assume that the potential vorticity flux is 
down-gradient. This assumption is a good one for barocli- 
nically unstable waves as seen from observations and theory 
(Edmon et al . , 1980) . For quasi-geostrophic motion this 
flux can then be written as [v*q*] = -K„ ** [q]/3y where K CT 
is a mixing coefficient of potential vorticity. Here we 4 
assume that Kq=KQsin y' so that no potential vorticity 
flux occurs at the sidewalls. This choice maximizes the 
mixing in the center of the channel where the meridional 
temperature gradient is strongest, i.e., where we expect 
maximum eddy activity. The amplitude. Kg, of the mixing 
coefficient is determined by substituting the down— gradient 
q-flux into the constraint ( 8 ) . The mixing coefficient 
depends on the strength of the heat flux and a volume- 
average potential vorticity gradient: 

fLy H 

(f/N ) 2 [v*3J^/3z]| dy 
K 0 = J 0 0 

f H fLy 

J J Mq]/ay Sin y» dy dz 

After finding Kg, which is dependent on the evolving 
fields, we can find [v* 5 *] from (7) . Treating these 
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quantities as time, as well as zona l . means, we can find an 
expression for & . T he f lux [v*$*] includes both the 
transient eddy flux, [v'S‘], and stationary eddy flux, 
[v*3*]. Of course, we are only parameterizing the former 
quantity since the latter is explicitly computed in the 
model as part of the predicted field of yJ . 

d. Condensational heating 

The flux of water vapor by transient eddies is v|m' 
where m is the specific humidity. The quantity, m, is a 
strong function of temperature and here we assume m=m(T) . 
That is, we assume that variations of specific humidity in 
a transient eddy_are entirely dependent on temperature so 
that m' = m(T)-m(T)=R h (m s (T)-m s (T) ) where R h is the relative 
humidity and m s is the saturated value of m. Then, m s can 
be found from the Clausius-Clapeyron equation. The relative 
humidity is assumed to depend on geographic location only 
so that we can introduce a global-scale variation repre- 
senting the continent-ocean moisture contrasts. 

The water vapor flux can be related to the eddy heat 
flux by expanding m s ' in terms of temperature deviation so 

that _ 

v'm' = R h Bm s (T) T'v' . 

&T 

A similar expression was derived by Leovy (1973) for use in 
planetary atmospheres modelling. 

We are interested in modelling latent heat release in 
regions of precipitat ion. Thus, we only consider regions 
where convergence of v'm' occurs, implying a loss of water 
vapor through condensation. In these regions, latent heat 
is released and we obtain only positive contributions to 
diabatic heating. We implicitly assume that evaporation 
and associated cooling has previously occurred at the ocean 
surface and does not affect the mid— latitude planetary- 
scale circulation. The contribution to diabatic heating is 

gLv 

q = (|V* v'm' | - ▼* v'm') 

2fC p T 

where C p is the specific heat of moist air and Ly latent 
heat of condensation. Note that Q c is zero or positive. 

7. Model equations 

Without further approximation, we would obtain a set of 
12 coupled nonlinear, first-order differential .equations 
for the evolution of the components of rg and r> p. However, 
these equations can be simplified by scaling (5) -(6) and 
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neglecting small terms. First, scale Y by UL and x and y 
by L where U' w 20 m sec“l and L^6000 km for these scales of 
motion. The nonlinear and local derivative terms in (5) 
are small (0(U//9L 2 )=0.03) compared to the advection of 
planetary vorticity by these waves (the term) . In 

(6) the terms involving the relative vorticity are also 
small by the same reasoning. However, the terms involving 
A must be retained since U A 2 /^5 is 0(1). 

Thus, we arrive at six predictive equations for the 
spectral components of and six diagnostic equations 
relating Jk b and it is not instructive to present the 
equations for the individual spectral components. However, 
it is relevant to consider the basic forms in terms of the 
zonal mean and wave components. 

The zonal-mean form of (5) is 

<^ > - kf d 2 [^ 0 ] - f [J(^ Q *,h)] = 0 

H dy 2 

where t P 0 = ^(z=0) . This equation expresses a balance 
amongst transient eddy vorticity flux, Ekman dissipation of 
the zonal mean flow and mountain torque induced by the 
planetary wave. The wave form of (5) is 

= “kf pYo* ” 1 J([^ 0 ]>h) 

“ ^x H H 

which represents a balance amongst planetary vorticity 
advection by the waves, Ekman dissipation and forcing by 
zonal flow over topography. 

The zonal-mean form of (6) is 
Atfe + [J(^B*^T*)] = -Jif d 2 [^ 0 ] Fi(0) + 

A l H dy 2 

_£ [J(^o*' h )3 + (f/NA ) 2 <[Q]dF;,/dz> 

A l H 

expressing changes in the zonal-mean temperature by plane- 
tary wave heat fluxes, adiabatic heating/cooling caused by 
zonal-mean vertical velocity that is induced by Ekman pump- 
ing and mountain torque, and zonal-mean diabatic heating. 
The contribution from the vertical variation of the eddy 
vorticity flux in (6) is small compared to the other terms 
and has been neglected. 

The wave form of (6) is 
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A.^t* + + J([^b 3*$r*) _ ^ A^jh* ~ 

St r Sx 

kf 7*^,* + f_J(t^c3/h ) + (f/N* ) 2 <Q*dF 1 /dz> 

A l H 

in which the change in waves' thermal field is determined 
by wave advection of the mean temperature field, mean flow 
advection of the wave temperature field, adiabatic heating/ 
cooling induced by vertical variation in planetary vortic- 
ity advection, Ekman pumping and zonal-mean flow over 
topography, and, lastly, longitudinally-varying diabatic 
heating. The latter term includes damping by thermal 
radiation and transient eddy sensible heat flux and forcing 
by longitudinal variation in latent heating. 

These equations retain many fundamental feedbacks 
described in the conceptual model. The complete set of 
equations without the scaling approximations is also of 
interest and the simplified model should be compared with 
the complete formulation. For example, the simplified set 
does not allow for baroclinic instability of the planetary- 
scale waves. We have implicitly assumed that this process 
is not vital in generating the stationary planetary waves. 

8. Current status and future plans 

The code for numerically integrating the equations has 
been developed and is currently in the testing stage. Each 
process is being included in the model on an individual 
basis to test the code and examine the impact of the vari- 
ous processes on the circulation. 

In the future, the model will be tested with modifi- 
cations to the eddy flux parameterizations. It will also be 
expanded to allow a less severe spectral truncation. A 
seasonal variation of solar heating could also be investi- 
gated. 
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CONCLUSIONS AND RECOMMENDATIONS 


This report describes a model for the interaction of 
planetary and synoptic-scale waves based on some basic 
observations and theory. The intention is to investigate 
the effects of transient eddy (synoptic-scale) transports 
on the planetary-scale waves. Much of the synoptic-scale 
activity occurs over the northern oceans which are areas of 
sparse ground-based observations. Our knowledge of cumula- 
tive precipitation over these regions is unreliable, 
although its effect on the planetary-scale circulation may 
be significant on a climatological and seasonal basis. 

Models like the one described in this report will be 
useful in determining the impact of synoptic-scale waves on 
the larger-scale circulation. This model will complement 
other models with better spatial resolution and more ela- 
borate physical parameterizations. The strengths of the 
low-order model are in its ease of interpretation and 
numerical efficiency in long-term or extensive sets of 
experiments . 

The eventual goal is to ascertain and understand the 
dominant physical processes which operate in the extratrop- 
ical, large-scale circulation. In doing so, we can iden- 
tify what atmospheric parameters (e.g., precipitation over 
the oceans) must be observed by remote sensing techniques 
and what accuracy and frequency in the measurements must be 
achieved to allow successful modelling of the circulation. 
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APPENDIX 


Quasi-geostrophic Dynamics 

The quasi-geostrophic approximation to the equations of 
motion is suitable for motions whose vertical component of 
relative vorticity, as determined by horizontal velocity in 
a reference frame rotating with the Earth, is small com- 
pared to the Earth's vorticity (twice the angular velocity, 
n ). We scale the relative velocity by U and the hori- 
zontal gradients by L. Geostrophic balance between the 
Coriolis torque and pressure gradient force is the lowest- 
order diagnostic balance in the horizontal momentum equa- 
tions if 

U 

Rossby Number « l. 

2flL 

At the next order we obtain a predictive equation for 
the vorticity of the geostrophic (lowest-order) horizontal 
velocity: 

d ( V X fr+ f) = f M + 4* (A.l) 

dt a z 

where is the geostrophic streamf unction, f planetary 
vorticity=2 A sin ft, ft latitude, W vertical velocity, z 
vertical co-ordinate, 9 x y relative vorticity, t time and^ 
a source/sink of vorticity. The total derivative involves 
advection by the horizontal geostrophic motion only. Eq. 
(A.l) states that changes in absolute or total vorticity 
are caused by vertical stretching of vortex tubes (a result 
of angular momentum conservation in a rotating system) and 
sources/sinks. 

The potential temperature, 0, (the temperature an air 
parcel would have if brought to the surface adiabatically) 
is governed by the thermodynamic equation: 

d /g0'\ + N 2 W = Q (A. 2) 

dtle 0 1 

where g is gravitational acceleration, N=(g dlne 0 /dz)V2 
the Brunt -Vaisala frequency, ©'the deviation from the 
global mean 0 O , Q diabatic sources/ sinks. For hydrostatic 
motion (an excellent approximation for these motions) the 
streamf unction and 0' are, related by the thermal-wind equ.: 
gfi' = f 3^ 

®o * z (A. 3) 
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